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We have performed ®^'^^Ga nuclear magnetic resonance (NMR) and nuclear quadrupole reso- 
nance (NQR) and muon spin rotation/resonance on the quasi two-dimensional antiferromagnet 
(AFM) NiGa2S4, in order to investigate its spin dynamics and magnetic state at low temperatures. 
Although there exists only one crystallographic site for Ga in NiGa2S4, we found two distinct Ga 
signals by NMR and NQR. The origin of the two Ga signals is not fully understood, but possibly due 
to stacking faults along the c axis which induce additional broad Ga NMR and NQR signals with 
different local symmetries. We found the novel spin freezing occurring at T{, at which the specific 
heat shows a maximum, from a clear divergent behavior of the nuclear spin-lattice relaxation rate 
l/Ti and nuclear spin-spin relaxation rate I/T2 measured by Ga-NQR as well as the muon spin 
relaxation rate A. The main sharp NQR peaks exhibit a stronger tendency of divergence, compared 
with the weak broader spectral peaks, indicating that the spin freezing is intrinsic in NiGa2S4. The 
behavior of these relaxation rates strongly suggests that the Ni spin fluctuations slow down towards 
Tf, and the temperature range of the divergence is anomalously wider than that in a conventional 
magnetic ordering. A broad structureless spectrum and multi-component Ti were observed below 
2 K, indicating that a static magnetic state with incommensurate magnetic correlations or inho- 
mogeneously distributed moments is realized at low temperatures. However, the wide temperature 
region between 2 K and Tf, where the NQR signal was not observed, suggests that the Ni spins 
do not freeze immediately below Tf, but keep fluctuating down to 2 K with the MHz frequency 
range. Below 0.5 K, all components of 1/Ti follow a behavior. We also found that l/Ti and 
I/T2 show the same temperature dependence above Tf but different temperature dependence below 
0.8 K. These results suggest that the spin dynamics is isotropic above Tf, which is characteristic of 
the Heisenberg spin system, and becomes anisotropic below 0.8 K. 



PACS numbers; Valid PACS appear here 

I. INTRODUCTION 

Recently, vigorous theoretical and experimental stud- 
ies have been performed on various compounds with geo- 
metrically frustrating lattices.^ Among such compounds, 
a two-dimensional (2D) triangular lattice is quite intrigu- 
ing because it is the simplest and most fundamental 
structure. Quite recently, a new quasi 2D triangular anti- 
ferromagnet NiGa2S4 was discovered as the first example 
of a bulk low-spin antiferromagnet with an exact trian- 
gular lattice.^ NiGa2S4 is a layered compound with the 
central NiS2 block layers of edge-sharing NiSe octahcdra 
and top and bottom sheets of GaS4 tetrahedra (Fig. 1). 
These slabs are stacked along the c axis and connected 
with each other by a weak van der Waals force. Since 
the NiS2 layers are effectively decoupled, and the Ni- 
Ni distance along the c axis is more than three times 
longer than the Ni-Ni distance along the a axis, NiGa2S4 
has been regarded as a nearly ideal 2-D triangular sys- 
tem. The electronic configuration of magnetic Ni^+ (3d*) 
ions is ^2g;Sg — !)■ Despite strong antiferromagnetic 
interactions (the Weiss temperature 0w — —80 K), it 
was reported that no long-range magnetic order is de- 



tected down to 0.35 K by susceptibility, specific-heat, 
and neutron-diffraction measurements.^ Instead, incom- 
mensurate short-range order with nanoscale correlation, 
whose wave vector is q = {rj, rj, 0) with = 0.158 ~ 1/6, 
was revealed by the neutron diffraction. In addition, a 
quadratic temperature dependence of the specific heat 
below 4 K indicates the presence of coherent gapless lin- 
early dispersive modes at low temperatures, although the 
bulk measurements show a small anomaly around 10 K 
{— Tf) Ref.[2]. However, the origin of the dependence 
of the specific heat (C) has not been understood yet. 

In this paper, we report the results of Ga nuclear- 
magnetic-resonance (NMR) and nuclear-quadrupole- 
resonance (NQR) measurements on NiGa2S4, which have 
been performed in order to investigate from a micro- 
scopic viewpoint the origin of the anomaly at Tf = 10 
K, spin dynamics at low temperatures, and the mag- 
netic ground state. We also report the results of muon- 
spin-rotation/relaxation (//SR) measurements, which de- 
tect the slow dynamics beyond the NMR limitation. We 
found that the nuclear spin-lattice relaxation rate 1/Ti, 
nuclear spin-spin relaxation rate 1 /r2 and the muon spin 
relaxation rate A all exhibit a clear divergence in ap- 



2 



preaching from above Tf ~ 10 K. In addition, a broad 
structureless Ga-NQR spectrum and inhomogeneous dis- 
tribution of Ti were observed below 2 K. These results 
strongly suggest that the Ni moments freeze out below 
Tf and give rise to inhomogeneous internal fields at the 
Ga sites below 2 K. Moreover, it was found that the 
spin dynamics is isotropic above Tf , which is shown from 
the identical temperature dependence of 1/Ti and I/T2. 
Isotropic spin dynamics is considered to be characteris- 
tics of Heisenberg spin system. The different temperature 
dependence was observed at low temperatures between 
1/Ti and I/T2: 1/Ti follows a T^ dependence below 0.5 
K, and I/T2 is almost linear in T below 2 K. The dif- 
ferent temperature dependence suggests anisotropic spin 
dynamics at low temperatures. 




FIG. 1: (Color online) Crystal structure of NiGa2S4. The 
lattice constants are a = 3.624 A, c = 11.996 A. 



II. EXPERIMENTAL 

Polycrystalline and single-crystal samples of NiGa2S4 
were synthesized as described in the literature.^ Pow- 
der x-ray measurements at room temperature and neu- 
tron diffraction measurements in the temperature range 
between 1.5 K and 300 K confirmed that NiGa2S4 re- 
tains the trigonal crystal structure down to 1.5 K with 
P3ml symmetry.^ Three samples (polycrystalline and as- 
grown single crystalline samples of nominal NiGa2S4, and 
a polycrystalline sample of NiGa2S4.o4 ) were used in 
the measurements. The same batch of polycrystalline 
NiGa2S4 used in the neutron diffraction^ was measured 
in the present NQR experiments, and the single crystal 
was used in NMR. Since NQR results (NQR spectrum 
and 1/Ti) are essentially the same in the three samples, 
we consider that the NQR results are determined by the 
intrinsic magnetic properties, which are not affected by 
a spurious impurity phase. 

The NMR/NQR measurements were performed on two 
Ga isotopes [^^^Ga (/ = 3/2): = 10.219 MHz/lOkOe 



and = 0.19 X lO-^^cm^, '^iGa (/ = 3/2) : ■^^7 = 
12.984 MHz/lOkOe and ^^Q = 0.16 x lO^^-^cm^, where 
7 and Q are the gyromagnetic ratio and the electric 
quadrupole moment, respectively]. NMR/NQR spectra, 
1/Ti and I/T2 were measured by Ga spin-echo signals in 
the temperature range of 75 mK - 200 K. Zero-field and 
longitudinal-field /xSR measurements in the temperature 
range of 1.8 K to 250 K were performed at the 7rA-port 
of the Meson Science Laboratory at KEK in Tsukuba, 
Japan. The powder sample was attached to a silver "cold 
plate" by GE-varnish. A ^He gas flow cryostat was used 
for the /iSR measurements. 



III. EXPERIMENTAL RESULTS 

A. NMR and NQR spectra 

NMR spectra of a single-crystalline NiGa2S4 were ob- 
tained by sweeping the external field. Figure 2 (a) dis- 
plays the Ga-NMR spectra against the field at 40, 17 and 
1.6 K. Here, the NMR frequency is fixed at 93.5 MHz 
and the external magnetic field is applied parallel to the 
c-axis. 

In general, NMR spectra for nuclear spin I{~ 3/2) 
show three peaks, which are composed of an intense cen- 
tral peak arising from the 1/2 <-!■ —1/2 transition and two 
satellite peaks from 3/2 1/2 and — 1/2 ^ —3/2 tran- 
sitions when electric field gradient exists at the observed 
nucleus site. Since NiGa2S4 has one crystallographic Ga 
site surrounded by four S atoms, which forms a GaS4 
tetrahedron, three peaks should be observed for each of 
the ^^Ga and ^^Ga nuclei. However, the spectrum at 40 
K exhibits asymmetric center peaks and two asymmetric 
satellite peaks, indicative of the existence of two Ga sites 
(Ga(l) and Ga(2) sites) with different EFGs. With de- 
creasing temperature, the intensity of the NMR signals 
from the Ga(l) site decreases, and NMR signals from 
only one Ga site were observed at 17 K as shown by the 
dark gray area in Fig. 2 (a). It should be noted that 
measurement of T2 on NMR spectra was needed in or- 
der to discuss the intrinsic intensity ratio of the two Ga 
sites. By further cooling, all NMR signals disappeared 
below Tf ~ 10 K, and extremely broad NMR spectrum 
reappears below 2 K as shown in the bottom of Fig. 2 
(a). The extremely broad Ga-NMR spectrum strongly 
suggests the occurrence of static magnetism with inho- 
mogeneous internal field well below Tf. 

The existence of two Ga sites with different EFGs, 
which was indicated by the Ga NMR, was confirmed by 
the NQR measurement. Ga-NQR spectra from polycrys- 
talline sample were obtained by measuring the spin-echo 
intensity as a function of frequency. Unlike the NMR, 
the NQR spectrum for the 7 = 3/2 nuclei should con- 
sist of a single peak of the ±1/2 ^ ±3/2 levels. We ob- 
served two intense peaks with narrow width and two weak 
peaks with broad width as shown in Fig. 2 (b). The fre- 
quency ratio of two intense peaks ( 16.5 MHz / 10 MHz) 
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FIG. 2: (a)'^^'"Ga-NMR spectra from single-crystal NiGa2S4 
at 40 K (curve in upper part), 17 K (gray area), and 1.6 K 
(lower part). NMR intensity normalization, i.e. taking into 
account difference in the T2 relaxation rate {T2 correction) 
has not been made. The resonant frequency is 93.5 MHz and 
the magnetic field is applied along the c-axis. The time be- 
tween 1st and 2nd pulses to observe spin echo is 30 /isec. The 
^'■^^^^ Cu signals are from an NMR coil, (b)""'^^ Ga-NQR spec- 
tra after T2 correction at 40 and 1.5 K using poly crystalline 
samples. NQR intensity is normalized by temperature. 



is nearly the same as that of two weak peaks (12.5 MHz 
/ 7.5 MHz), which is equal to the ratio of the nuclear 
quadrupole moment {^^Qf^^Q = 1.589). This indicates 
that there are two sets of distinct NQR signals, and re- 
veals the existence of two Ga sites with different EFGs 
in NiGa2S4. From the separation between two satellite 
peaks in the NMR spectra at 40 K, we deduced the NQR 
frequency along the principal axis (c axis) of the EFG at 
the Ga site {i^z)- We evaluated the asymmetry parameter 
V ~ i^v ~ ^x)lvz using the relation 



tion is ~ 3.96, suggesting that the S deficiency is at most 
~ 1%, and that the configuration of the triangular struc- 
ture is rather good. At present, the origin of the two Ga 
sites is not identified, but we point out possible inclusion 
of different stacking units, closely related to the struc- 
ture of NiGa2S4, might give rise to an additional Ga site 
with a different EFG. A tiny amount of sulfur deficiency 
and/or disorder is considered to exist in the 'outer' sulfur 
layer, which is shown by S(l) in Fig. 1. It is considered 
that the sulfur deficiency and/or disorder is more eas- 
ily introduced in the outer S(l) layer than in the inner 
S(2) layers which are coupled strongly with Ni^+ ions 
and form the NiS2 block layer. The fraction of the two 
Ga sites estimated from the Ga-NQR intensity at 40 K, 
which is normalized by I/T2 values, is Ga(l) : Ga(2) = 
0.78 : 0.22. If we assume that the Ga(2) NQR is ascribed 
to the Ga site which is influenced by the sulfur deficiency 
at the S(l) site, we estimate that the S(l) deficiency is ap- 
proximately 7 % from the intensity ratio of two Ga NQR 
signal because one S(l) deficiency gives an influence to 
three Ga atoms in the low concentration limit. However, 
this possibility might be excluded from the ICP result. 
Instead, we consider that the stacking faults along the 
c axis, which result in a different stacking unit from the 
bulk NiGa2S4, might be the origin of the Ga(2) site. In 
any case, the linewidth of the Ga(2) NQR signal is twice 
broader than that of the Ga(l) signal, it is reasonably 
considered that the Ga(l) signal arises from the Ga site 
with the regular crystal structure, and the Ga(2) signal 
arises from the Ga site with disorder and/or defects in 
the structure. This assignment is consistent with the fact 
that r\ at the Ga(2) (~ 0.75) is larger than r\ at the Ga(l) 
(- 0.37). 

The NQR-signal intensity at both Ga sites decreases 
with decreasing temperature and the NQR signals disap- 
pear around Tf ~ 10 K, indicative of a magnetic anomaly. 
On further cooling, enormously broad NQR spectra were 
observed as in the NMR spectra below 2 K. Since no obvi- 
ous structure was found in the spectra, it is inferred that 
the internal field at the Ga sites is widely distributed in 
a static magnetic state at low temperatures. Such mag- 
netic state is discussed in section IV B, C. 
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where vq is the NQR resonance frequency in Fig. 2 (b). 
The values of :/q, v^, and 77 are listed in Table I. 

It is noted that two Ga NQR signals are observed, 
even though there exists only one crystallographic site 
for Ga in perfect NiGa2S4. From careful Inductively 
Coupled Plasma (ICP) and scanning electron microscope 
(SEM) measurements^ it was revealed that S occupa- 



TABLE I: The quadrupole parameters for ''^'^^Ga(l) 
^®'^^Ga(2) derived from the peak position of the NMR and 
NQR spectra. The ratio of quadrupole moment of Ga iso- 
topes which corresponds to that of NQR resonance frequency 
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B. Knight Shift 

Next, we show the temperature dependence of the 
Knight shift (K) measured at the central peaks of the 
Ga-NMR spectrum. In the K estimation, the shift orig- 
inating from the 2nd-order quadrupole interaction was 
subtracted. Figure 3 displays the temperature depen- 
dence of K{T) for 7iGa(l) and ^^Ga(2), along with the 
behavior of the bulk susceptibility normalized by the be- 
havior of the Knight shift above 100 K.^ In general, K{T) 
is decomposed as. 
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where Kspin(T) and Ko^h are the spin and orbital parts 
of Knight shift, respectively. -fCgpin is related to the bulk 
susceptibility xC^) originating from the Ni spins as 



FIG. 3: Temperature dependence of the Knight Shift K for 
■^^Ga(l)(») and ■^^Ga(2)(o) sites. The dotted curves in the 
main panel give the temperature dependence of susceptibility 
X normalized to K above 100 K. Inset: K vs. x pld for the 
two ^'^Ga sites. Dotted lines are the linear fits between 70 K 
and 200 K. 



where Am, -/Va, and /j,b are the hyperfine coupling con- 
stant between the Ga-nuclear spin and Ni-Sd spins, Avo- 
gadro's number and the Bohr magneton, respectively, 
i^orb is related with the Van Vleck susceptibility, which 
is temperature independent in general. 

As shown in Fig. 3, the Knight shift at both sites 
follows the bulk susceptibility above 80 K. The hyper- 
fine coupling constant at both sites [17.70 ± 0.20 
(kOe//iB) for Ga(l), 7.67 ± 0.36 (kOe//iB) for Ga(2)] 
was evaluated from the slopes in the K-x plot displayed 
in the inset of Fig. 3. However, K at the Ga(l) site 
could not be measured below 50 K due to the decrease 
of the NMR-signal intensity, and K at the Ga(2) site 
deviates from the temperature dependence of x(T) be- 
low 80 K. The Knight shift probing the microscopic 
susceptibility decreases although the bulk susceptibil- 
ity continues to increase. It seems that a precursor 
of the magnetic anomaly starts below 80 K. This cor- 
responds well to the Weiss temperature obtained from 
x(T). We point out that a similar deviation of K{T) 
from x(T) at low temperatures was observed in a num- 
ber of geometrically frustrated spin systems such as 
SrCr8Ga40i9(SCGO),4 Ba2Sn2ZnCr7pGaio-7p022^ and 
FeSc2S4.^ In these compounds, the possibility of spin sin- 
glet state and/or gapped state with the defect like contri- 
bution has been suggested. However such possibilities are 
not applicable to the present case because the divergence 
of 1/Ti, which is discussed later, suggests a magnetic 
ground state. In the magnetic state, there exist static 
moments giving rise to inhomogeneous internal fields at 
the Ga sites. The different temperature dependence be- 
tween microscopic and macroscopic spin susceptibilities 
in the highly frustrated compounds, which are probed 
with NMR and bulk susceptibility respectively, remains 
unclear, and deserves to be understood theoretically. 



C. Nuclear spin-lattice and spin-spin relaxation 
rate, l/Ti and I/T2 

The nuclear spin-lattice relaxation rate (longitudi- 
nal rate) l/Ti was measured by the saturation-recovery 
method with saturation pulses. Above Tf, the recov- 
ery curve m{t) of the nuclear magnetization M{t) at 
time t after the saturation pulses, which is defined by 
m{t) = [M{oo) - M(t)]/M((X)), is consistently fitted 
by the m{t) oc exp(— 3t/Ti) at both Ga sites as shown 
on a semilogarithmic scale in Fig. 4 (a). Thus a sin- 
gle Ti component was determined above 10 K. The nu- 
clear spin-spin relaxation rate (transverse rate) 1 IT2 was 
measured by recording the spin-echo intensity I{2t) by 
changing the time separation r between the 7r/2 (1st) and 
TT (2nd) pulses. /(2r) was well fitted by the relation of 
/(2r) oc exp(— 2r/T2) at both Ga sites as shown in Fig. 4 
(b). Here, the pulse length of a 7r/2 is approximately 5 
/isec. 

Between T[ and 2 K the NQR signals are not observed 
due to the extremely short Ti and T2. This is due to 
the critical slowing down of the spin fluctuations to the 
NQR frequency scale, which results in the rapid spin re- 
laxation. On further cooling, a broad NQR spectra in 
Fig. 2 (b) was observed below 2 K. The m{t) exhibits 
the multi-component behavior, which is characterized by 
the upward curve as shown in Fig. 5 (a), which m(t) is 

consistently fitted by the relation m{t) oc exp{~^3t/Ti) 
in Fig. 5 (a). This is shown by the straight line in the 
semi-log plot between m{t) and ^/t. This relation is often 
observed when 1/Ti is inhomogeneously distributed due 
to the presence of the relaxation center, e.g. a amount of 
doped magnetic impurities and/or crystal imperfectionii 
In such a case, 1/Ti is normally adopted as 1/Ti. On 
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the other hand, it was found that /(2t) follows the same 
7(2t) oc exp(— 2r/T2) relation even below Tf. The in- 
set of Fig. 6 shows the plot of the /(2r) against 2r at 
76, 220, and 485 mK. Thus single-component I/T2 was 
determined down to the lowest temperature. 

The temperature dependences of 1/Ti, 1/Ti and I/T2 
determined by above fitting procedures are displayed in 
Figs. 7 (a) and 7 (b), together with the temperature 
dependence of the integrated NQR intensity multiplied 
by temperature, which is normalized to the value at 100 
K [Fig. 7 (c)]. A T2 correction has been made for the 
estimation of the NQR intensity in Fig. 7 (c), and we 
assume that T2 is independent of frequency below 2.5 K. 
Figures 7 (a) and (b) display the temperature dependence 
of 1/Ti and l/r2 measured at the two ^^Ga sites above 
10 K, 1/Ti measured at various frequencies, and I/T2 
at 12 MHz, respectively. The ratio of 1/Ti between the 
two isotopes at high temperatures is (l/^^Ti)/(l/^^Ti) 
= 0.67 ± 0.05, which is in good agreement with the ratio 
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distributed as shown in Fig. 5 (b). (c) Temperature depen- 
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of gyromagnetic ratio (^^7/^^7)^ = 0.62. This indicates 
that 1/Ti is dominated by the magnetic interaction with 
the electrons, not by the electric quadrupole interaction. 
Above 80 K, l/Ti of ^^Ga(l) is nearly constant, which 
is often observed in a localized-moment system when the 
temperature is higher than the exchange energy between 
the localized moments. Below 80 K, as shown in Fig. 7 
(a), 1/Ti gradually increases with decreasing tempera- 
ture, and correspondingly the intensity of the NQR signal 
decreases. Such behavior is considered as a precursor of 
the spin freezing of Ni-Sd spins. If we compare the tem- 
perature dependence of 1/Ti at the two Ga sites, ^^Ga(l) 
and ^^Ga(2) sites, the magnitude of 1/Ti at the Ga(l) 
site is three times larger than 1/Ti at the Ga(2) site. It 



should also be noted that the onset temperature of the 
divergent behavior below about 80 K is higher at the 
Ga(l) site than at the Ga(2) site. From the comparison 
of the temperature dependence of 1/Ti at the two Ga 
sites, we conclude that the magnetic anomaly seen in Tf 
is an intrinsic nature of NiGa2S4, because the Ga(l) site 
possesses a narrower NQR spectrum, indicative of the ho- 
mogeneous site. If Tf were induced by inhomogeneity of 
the sample, such as sulfur disorder and/or stacking faults, 
l/Ti at the Ga(2) site, where larger inhomogeneity is sug- 
gested by the NQR spectrum, would show the magnetic 
precursor from higher temperatures. Obviously, this is 
not the case. A similar difference between the Ga(l) and 
Ga(2) sites, which indicates that the magnetic anomaly 
at Tf is intrinsic in NiGa2S4, was also observed in I/T2 
and NQR-intensity results as shown in Figs. 7 (b) and 
(c). 
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As noted above, in the temperature range between 
10 and 2 K, the NQR signal was not observed due to 
the extremely short Ti and T2 beyond the limit of the 
NMR measurement. Below 2 K, the relaxation rate is 
distributed inhomogeneously in the low-temperature re- 
gion. This multi-component Ti behavior is often ob- 
served in the spin-glass® and Kondo-disordered systemsi^ 



1/Ti determined from m(i) cx exp ^—y3t/Tij decreases 

strongly with decreasing temperature as shown in Fig. 7 
(a). Here, Ti/3 is the time when m{t) decays to 1/e, is re- 
garded as a characteristic value of Ti in multi-component 
Ti. 1/Ti was measured at several NQR frequencies from 
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9.5 and 28 MHz at 1.5 K, but the values of Ti do not 
change so much in this frequency range. The tempera- 
ture dependence of 1/Ti below 2 K is shown in Fig. 7 (a) 
in the frequency range between 9.3 and 13 MHz, and is 
independent of the frequency. 1/Ti measured in 10 kOe 
shows the same temperature dependence as that in zero 
field, although the value of 1/Ti in 10 kOe is approxi- 
mately six times smaller than l/Ti in zero field. This is 
due to the difference of the matrix elements of the tran- 
sition probability. 

We found that 1/Ti is proportional to between 0.8 
and 0.1 K. It should be noted that, although the spin 
dynamics is highly inhomogeneous, the dependence 
holds in all 1/Ti components, because the experimental 
data of m{t) at several temperatures below 0.5 K lie on 
the same curve when plotted against tT^ as shown in 
Fig. 5 (b). Therefore, the dependence of l/Ti is in- 
trinsic behavior at low temperatures in NiGa2S4. Tem- 
perature dependence of 1/Ti is discussed in Sec. IV. 

Figure 7 (b) shows temperature dependence of I/T2. 
I/T2 at both Ga sites also diverges at T{, but I/T2 at 
the Ga(l) site has larger values and starts to diverge at 
a higher temperature, which is in agreement with the 
l/Ti results. Below 2.5 K, I/T2 follows a T-linear tem- 
perature dependence down to the lowest temperature as 
shown from the plot of /(2t) against 2rT at several tem- 
peratures in the main panel of Fig. 6. The large I/T2 val- 
ues and the strong temperature dependence of 1 /T2 that 
follows 1/Ti behavior above T{ indicate that I/T2 is not 
determined by the nuclear dipole interaction, but by the 
Ti process. We found that the ratio of T2/T1 in NiGa2S4 
is 4.2 ± 0.2 above Tf . We show temperature dependence 
of 1/Ti and I/T2 divided by 4 in Fig. 8. It should be 
noted that l/(4r2) follows quantitatively the same tem- 
perature dependence as 1/Ti above Tf for both Ga sites, 
and that different temperature dependence is observed in 
1/fi and l/(4r2) below 0.8 K. In general, 1/Ti (I/T2) 
measured with the NQR spectra are affected by the spin 
dynamics perpendicular (parallel) to the principal axis 
of EFG. Therefore, it is considered that spin dynamics 
is isotropic above Tf, but becomes anisotropic below 0.8 
K. We suggest that the isotropic spin dynamics, charac- 
terized by the Heisenberg spin system, are slowing down 
and changed to the anisotropic ones below 0.8 K, where 
the in-plane spin dynamics is suppressed by the static 
magnetism, but the out-of-plane spin dynamics remains 
active and is proportional to temperature. Detailed tem- 
perature dependences above and below T{ are discussed 
in Sec. IV. 



D. Muon spin rotation/relaxation 

Muon spin rotation/relaxation (/iSR) experiments are 
particularly suitable for detecting slow spin dynamics 
with an extremely high relaxation rates beyond the NMR 
experimental limitation. Therefore, NMR and /iSR are 
complementary to each other. Detailed /zSR measure- 



ments in the temperature range between 2 K and Tf 
published in the literature.^" Here we concentrate on the 
zero -field /iSR (ZF-/iSR) experiments carried out at the 
Meson Science Laboratory KEK. 



NiGajS^ ZF 7- dependence 




t ( ^ sec ) 

FIG. 9: Zero-field ^SR asymmetry relaxation functions G{t) 
in NiGa2S4 at various temperatures. The dotted curves are 
fits to the power exponential form G{t) — ylexp [—(At)''] 

Figure 9 shows the ZF /iSR asymmetry relaxation func- 
tion G(t) in NiGa2S4 at various temperatures. The re- 
laxation function is gradually changing from Gaussian to 
Lorentzian by approaching Tf. The experimental asym- 
metry data were fitted to the "power exponential" form 

G{t)=Aexp [-(A tf], 

where A is the initial muon decay asymmetry, A is a gen- 
eralized relaxation rate, and the exponent (3 interpolates 
between exponential {(3 = 1) and Gaussian (/3 = 2) lim- 
its. This form is suitable for characterizing spin dynam- 
ics, i.e. Gaussian (exponential) function suggests that 
muons see static (dynamics) magnetic fields during their 
lifetimes (~ 2.2//s). /3 > 1 indicates an intermediate 
dynamics between two limits, and (3 < 1 suggests inho- 
mogeneous dynamics with multi-components. 

The relaxation rate A, exponent /?, and asymmetry A 
are plotted against temperature in Fig. 10. Above 100 
K, A is small and f3 is close to 2. Both indicate that the 
nuclear dipole effect is dominant in the relaxation. With 
decreasing temperature, f3 gradually decreases from 2 to 
1, and A increases, although A is nearly constant down to 
Tf. These results indicate that fiSK can detect the slow 
spin dynamics which is beyond the NMR limitation, since 
the NMR intensity gradually decreases in this tempera- 
ture range. A shows a distinct peak at Tf ~ 10 K, below 
which the asymmetry A decreases to approximately 1/3 
of the values above Tf. As evident from Fig. 9, the so- 
called "1/3 tail" is lifted up at 2 K, which indicates that 
the static field along the implanted muon direction ap- 
pears at least 2 K. In general, the muon spin polarization 
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FIG. 10: Temperature dependence of (a) /iSR relaxation rate 
(A), (b) exponent /3 and (c) asymmetry A derived from the 
fitting of zero-field asymmetry by G{t) = Aexp [—(A t)^]. 



IV. DISCUSSION 
A. Spin dynamics above Tf 

In this section, we discuss spin dynamics in NiGa2S4 
in the paramagnetic state, which are revealed by l/Ti 
and I/T2 at the Ga(l) site, because the Ga(l) is consid- 
ered to be more homogeneous than the Ga(2) site from 
the narrower NQR spectrum shown in Fig. 2 (b) . Above 
Tf, 1/Ti and 1/(4T2) behave quantitatively the same as 
shown in Fig. 8 and 11, and A derived from the /xSR ex- 
periments also follows the same temperature dependence. 
At temperatures higher than 80 K, 1/Ti and 1/(4X2) are 
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function below Tf is characterized by rapid depolarization 
of 2/3 of the initial polarization, followed by slow dy- 
namic relaxation of the remaining 1/3 component. This 
is a characteristic signature of a highly disordered mag- 
netic state in which the moments are quasi-static on the 
time scale of the muon lifetime 1^ In the present mea- 
surement, the rapid depolarization could not be detected 
due to the high inhomogeneity of the static moments, but 
the 1/3 component was detectable. It is considered from 
the residual asymmetry below Tf that almost the entire 
region in the sample is in the magnetic state. 

It should be noted that relaxation of the 1/3 tail was 
observed at T = 2 K. Therefore, strong relaxation demon- 
strates the existence of slow dynamics even below Tf, 
which is consistent with the experimental fact that the 
Ga NQR spectrum is not observed down to 2 K. The 
present /iSR measurements suggest that magnetic fluc- 
tuations of Ni spins are gradually slowing down below 
80 K, and that the inhomogeneous static field appears at 
the implanted muon site at low temperatures. However, 
fluctuations of the Ni spins remain strong even far below 
Tf . Spin dynamics is discussed in the next section on the 
basis of NMR and /LtSR relaxation behavior. 



FIG. 11: Temperature dependence of l/Ti and l/(4r2) of the 
Ga(l) NQR spectra, and /iSR relaxation rate A. The inset 
shows the plot of l/Ti and l/(4r2) against (T - rf)/rf. 



nearly constant due to local-moment fluctuations by the 
exchange interaction between the Ni spins. Below 80 K, 
1/Ti and 1/(4T2) start to diverge towards Tf. As seen in 
the inset of Fig. 11, the divergent behavior of l/Tf and 
1/(4T2) is approximately fitted by 

1/Ti,2 cx e-" 

with £ = (T — Tf)/Tf, and Tf and i' are estimated to be 
Tf = 10.3±1.0, andt/ = 0.45±0.10. According to theoret- 
ical prediction, i' is 0.3 and 0.8 in 3-D and 2-D Heisenberg 



systems, respectively 
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0.5 suggests that magnetic 



anomaly might occur due to weak interactions along the 
c axis, since magnetic order is not realized at a finite 
temperature in the 2-D Heisenberg system. It should be 
noted that the critical divergence of the relaxation rate 
was observed over an unusually wide temperature range 
0.2 < e < 3, in which the 2-D short-range correlation is 
highly developed. This is in marked contrast with the or- 
dinary critical behavior in the 2-D Heisenberg systemic, 
for which the critical divergence is limited to the vicinity 
of T/v (e < 0.5). Such a remarkable divergent behavior 
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persisting over a wide temperature range seems to be one 
of the key features of the 2-D triangular lattice, because 
similar critical behavior over wide temperature ranges 
was reported in LiCr02 and HCr02 with 2-D triangular 
lattices J^iii 

In addition, it was found that 1/{AT2) shows quan- 
titatively the same temperature dependence as 1/Ti. 
This indicates the spin dynamics is isotropic. When the 
nuclear dipole interaction is negligibly small, there are 
the relation between l/Ti and I/T2 measured with the 
NQR signal arising from the ±1/2 ^ ±3/2 of / = 3/2 
transitions 



1 
1 



Here, Ga{ui) is the spectral density of the longitudinal 
(a = z) and transverse (a =±) components of the fluc- 
tuating local field^, and is given by: 



21 
2 



{SHa{t)SHa{0)) exp {iujt)dt, 



where 7 is the gyromagnetic ratio, and 6Hz{t) {dH±{t)) 
is the longitudinal (transverse) component of the fluctu- 
ating local fleld. If we assume an exponential correlation 
function, 



using the correlation time Tc of the fluctuation, 1 /Ti and 
I/T2 are expressed as follows. 



1 



1 + (ojTcy 



As seen in the equations, 1/Ti and I/T2 are determined 
by the dimensionless value oi ujTc, and are proportional 
to Tc (l/w^Tc) in the case of ujTc <^ 1 {ujTc 1). When 
Tc becomes longer due to the slowing down of the Ni 
spins with approaching Tf, the relaxation rates once in- 
crease with a maximum at luTc — 1 and then decrease 
with decreasing temperature. The experimental facts 
that I/T2 ~ 4(1/Ti) and both show quantitatively the 
same temperature dependence each other imply that the 
magnetic fluctuations are in the regime luTc <ti 1 and 
isotropic, and that the relation between 6H± and SHz is 
{{SHzY') ~ l.b{{5H±Y) in the temperature region above 
Tf. The relaxations behaviors are consistent with the 
Heisenberg spin system. 

With decreasing temperature, magnetic fluctuations 
are slowing down below 80 K, and seem to be static be- 
low Tf. Now, we discuss the temperature dependence of 
the correlation time Tc(T) of the magnetic fluctuations 



quantitatively using the experimental data for 1/Ti and 
ylhf at the Ga(l) site and magnetic susceptibility x- In 
the case when the nuclear spin-lattice relaxation rate at 
the Ga site is dominated by magnetic fluctuations of the 
Ni localized spins, 1/Ti is expressed as follows; 



1 



7. 



2X"{q.oj) 



where z{= 3) is the number of the nearest neighbor Ni 
sites to a Ga site, A{q) is the g-dependent hyperfine cou- 
pling constant at a Ga site from a Ni spin, and ^"(g, w) 
is the dynamical susceptibility and the sum is over the 
Brillouin zone. At temperatures much higher than Tf, 
the spin dynamics is determined by independent Ni mo- 
ments, and the local-moment {q = 0) susceptibility is 
given by 



Xo{T) 



1 - iujTc{T) '' 



where xo is the magnetic susceptibility per Ni atom (emu 
/ Ni-atom). 

We consider that A{q) and the dynamical susceptibil- 
ity are isotropic because 1/Ti and 1/(4T2) show nearly 
the same behavior in this temperature region, and take 
A{q) ~ Am/z = 17.7/3 = 5.9 (kOe/Ms). Then 1/Ti is 
described in the regime of coTc <^ 1 as. 



1 __o llkBT 
Tl Ml 



Alao{T)TciT). 



Therefore, the characteristic energy of the spin fluctua- 
tions r/fce, which corresponds to h/[Tc(T)k-B], is given 
by 

T/kB = h/[Tc{T)kB] - 3 ^^"^f^'^ xo(r)Tir. 

Mb 

= 1.4 X 10^ TiTxoiVA (K). 



Figure 12 shows the temperature dependence of 
r(r)/fcB = h/[Tc{T)kB\. T(T)/kB is constant above 100 
K with the value of T{T)/kB ^ 20 K. This is in good 
agreement with the exchange interaction J between Ni 
spins, which is estimated to be J ~ 20 K from the rela- 
tion 



ik^Ow 
z'S{S + l)' 



where we use the Weiss temperature 0w of —80 K, and 
the number of the nearest neighbor Ni sites for a Ni ion 
z' of 6. 

It was found that r{T)/kB shows thermally activated 
decrease below 100 K, described by the relation of 

r(T)/fcB = Aexp(-t//r). 

From the inset of Fig. 12, A and U are estimated to be 
25 K and 40 K, respectively. Here, U corresponds to a 
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feasible that Tf is regarded as Tkm in NiGa2S4. 

We point out that the spin dynamics in NiGa2S4 pos- 
sesses the novel 2-D spin character, discussed in the 2-D 
Heisenberg compounds for nearly two decades As dis- 
cussed above, it was revealed that some kind of magnetic 
correlations start to develop below |0w| ^ 80 K, which 
is far above Tf, and continue to grow down to Tf. The 
development of magnetic correlations over wide temper- 
ature range and the occurrence of a magnetic anomaly 
well below 0w seems to be characteristic of frustrated 
magnetism, but it is still unclear what kind of magnetic 
correlations are developing below 9w- The magnetic cor- 
relations developing in the frustrated systems is an in- 
teresting issue to be studied from experimental and the- 
oretical points of view. 



B. Ga-NQR Spectrum below 2 K 



FIG. 12: Temperature dependence of the characteristic en- 
ergy of the spin fluctuations r{T)/kB (see in text). The inset 
shows the plot of r{T)/kB against 1/T 



binding energy of the local spin configuration determined 
by J. In general, we expect from the temperature de- 
pendence of 1/Ti that magnetic correlation length £,{T) 
develops exponentially with decreasing temperature, be- 
cause r(T) is related to ^(T) as r(r) cx [c/^(r)]" in 
the Heisenberg spin system. However, neutron exper- 
iments showed that the C(T) remains a couple of lat- 
tice constants even at 1.5 It seems that r(r) is 
not directly related to ^(T) in NiGa2S4. Alternatively, 
we point out that the similar thermally-activated be- 
havior was reported in ESR studies on 2-D triangu- 
lar compounds of LiCr02 and HCr02r^ and in NMR 
studies on a Heisenberg kagome lattice antiferromagnet 
KFe3(OH)6(S04)2'^ Particularly, Ajiro et al^ pointed 
out that the activation energy U in LiCr02 and HCr02 
takes a universal value equivalent toU = (2. 7±0. 3) JS*^ — 
iksTKM, where Tkm is a critical temperature for a 
Kosterlitz-Thouless-type phase transition. In the 2-D 
triangular spin system, the possibility of the "Z2 vor- 
tex state" has been suggestedJ^ Vortices related to the 
topologically stable defects appear at higher tempera- 
tures and starts to coupled each other below Tkm- If 
we apply the above analyses to NiGa2S4, J7 ~ 40 K 
corresponds to U ^ 2JS^ = 4fcBTf. This relation is 
quite different with the that observed in dilute-alloy spin 
glassesJ^ In these compounds, the correlation time Tc 
obeys the Arrhenius law Tc = tq exp(i?a/fcBT) with a 
very large activation energy Ea ^ 20fcBTg. This indi- 
cates that the freezing behavior is observed near Tg and 
the spins slow down immediately. We point out that the 
freezing behavior observed in NiGa2S4 is different from 
the canonical spin-glass systems, but the possibility of 
the "Z2 vortex" transition occurring in NiGa2S4. It is 



As shown in Fig. 2 (b), a broad and structureless spec- 
trum was observed at 1.5 K. This spectrum indicates the 
presence of inhomogeneous static magnetic fields at the 
Ga sites. When the internal magnetic field appears at an 
Ga nuclear sites below Tf with the electric quadrupole in- 
teraction, the Zeeman interaction from the internal fields 
is added to the total nuclear Hamiltonian. The Zeeman 
interaction is expressed as 

Ti-Z — "Infll ■ Hint, 

where 7n is the Ga nuclear gyromagnetic ratio and H-mt 
is the internal field at the Ga nuclear site. The struc- 
tureless broad spectrum is approximately reproduced by 
using the inhomogeneous internal field at the Ga nuclear 
site, as shown in the inset of Fig. 13. The average and 
distribution width of the internal field at the Ga site are 
approximately 0.5 T and 0.5 T, respectively, and the cal- 
culated Ga-NQR spectrum is shown by the dotted line 
in the bottom figure of Fig. 13. In the calculation, we 
assume that the inhomogeneous internal field has only a 
c-axis component for simplicity. Hi^t arises mainly from 
the transferred hyperfine field and the dipolar field from 
the Ni ions, which is expressed as 

where rrii is the magnetic moment at the i-th Ni site and 
Vi is the vector connecting the Ga site to the i-th Ni site. 
The average value of the ordered moments is estimated 
to be ~ 1.0/iB from the first term of the right-hand side 
in the equation above. It is difficult to estimate accu- 
rately the ordered moments and the magnetic structure 
from such a broad spectrum, since the magnitude of the 
dipolar fields highly depends strongly on the orientations 
of the moments. 
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FIG. 13: '^^'^^ Ga-NQR spectra at 40 K (upper panel) and 
at 1.5 K (bottom panel). The broad spectrum in the bot- 
tom figure is approximately reproduced by the presence of 
inhomogeneous internal fields pointing along the c axis. The 
distribution of the internal fields is shown in the inset of the 
bottom figure. 



C. Spin dynamics below Tf 

Since an extremely broad signal was observed below 
2 K, it is concluded that the static magnetic fields are 
present below 2 K. However, over the wide tempera- 
ture range below Tf, slow spin dynamics remains as 
probed by NQR and /iSR. Taking into account the multi- 
exponential behavior in the relaxation curves at low tem- 
peratures, it is concluded that a 2-D inhomogeneous mag- 
netic state is realized in this compound. 

We found the spin dynamics becomes anisotropic be- 
low 0.8 K from the different temperature dependence of 
l/Ti and l/iiT^) as shown in Fig. 8. 1/Ti follows a 
dependence in the temperature range between 0.5 and 0.1 
K, and this dependence holds for all the spin dynam- 
ics detected by the Ga NQR as discussed above. On the 
other hand, 1/(4T2) shows a T-linear dependence below 
1 K. If we take into account that 1/Ti and 1/(4T2) probe 
spin dynamics along different directions, we can conclude 
that G±{uo) and Gz{Q) show the and T- linear depen- 
dence, respectively. It is noteworthy that isotropic spin 
dynamics above Tf change to anisotropic ones below 0.8 
K. 

In the following, we discuss the temperature depen- 
dence of 1/Ti below Tf on the basis of a conventional 
two-magnon process in a triangular Heisenberg antifer- 
romagnet. 1/Ti determined by the two-magnon process 
was discussed by Moriyaf^ and was extended to the tri- 



angular lattice by Maegawa.''^^ 1/Ti by this process is 
expressed as 



1 TT r^m f 



n 






JuJO 







N{ujfduj, 



where 7e is the electronic gyromagnetic ratio, G^j is a 
geometrical factor. is the maximum frequency of the 
spin wave, ojq is the spin-anisotropy energy related to 
the internal field at the Ni sites, and N{ijj) is the state 
density of magnons. If we assume the long- wave length 
approximation, N{lo) in the 3-D and 2-D spin waves is 



N{uj) 



(w2 



7^ 



3-dimensional 
2-dimensionaL 



where lJcx is the exchange frequency between the Ni 
spins. Using N{uj), the temperature dependence of 1/Ti 
determined by the 3-D and 2-D spin waves is given by 



1 



cx T 5 

^l/3D JTo/T 



Trn 

T 



(e^ - 1)- 



and 



ocT^ / 

2D Jto/T 



T 



(e^ - 1)' 



respectively. Here, Tm = hwm/kB and Tq = hw^/k^. The 
temperature dependences of (1/Ti)3d and (1/Ti)2d are 
calculated using T^ - |6iw| = 80 K, and Tq = 0.4 K, 
which are shown in Fig. 14 by dotted and solid curves, 
respectively. 

1/Ti data in the temperature range below 1 K is 
roughly reproduced by both calculations using the gap 
magnitude of Tq = 0.4 K, however a deviation between 
the calculation and 1/Ti data is observed below 0.15 K. 
The deviation is considered to be due to the distribu- 
tion of To, which is suggested from the inhomogeneous 
internal field as discussed in Sec. IV B. It seems that 
the 2-D model is more suitable for interpreting the whole 
temperature dependence below Tf , because the deviation 
between the data and the 3-D calculation above 1 K is 
more significant than that in the 2-D one. 1/Ti from the 
3-D model is approximately a T^ dependence above 1 K, 
which is incompatible with the experimental result. The 
2-D spin-wave model seems to be consistent with the T^ 
dependence observed in specific-heat experiments above 
0.35 K,— which suggests the existence of linearly disper- 
sive modes. We point out that the 2-D spin-wave model 
also give the w-linear density of states. 
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FIG. 14: Temperature dependence of l/Ti calculated by the 
3-D and 2-D spin-wave models are compared with the exper- 
imental data (see text). 



The temperature dependence of 1/{4T2) is consistently 
interpreted by the dominant T-linear Gz (0) and the small 
dependence from the G±{uj) contribution discussed 
above. 1/{AT2) below 1 K is approximately expressed by 

which is shown by a dotted curve in Fig. 14. 

In the Heisenberg triangular magnet with the 120° spin 
structure, there are three kinds of low-energy excitation 
modes, viz. swinging fluctuations in the xy plane, rota- 
tional fluctuations around one of the moment direction, 
and vertical deviation from the xy plane (see Fig. 4 in 
RefHi). It is considered that the lowest excitation mode 
is the swinging fluctuation, because the fluctuation con- 
serves energyfi^ Taking into account that NiGa2S4 has an 
incommensurate short-range correlation with q close to 
(1/6, 1/6, 0), which is a 60° spin structure,^ it is reason- 
able to consider that NiGa2S4 has magnetic fluctuations 
similar to those in the triangular compounds. It should 
be noted that the out-of plane dynamics is homogeneous 
although the in-plane dynamics is inhomogeneously dis- 
tributed. It is considered that the out-of-plane dynamics 
is related to the lowest excitation, which is the swinging 
fluctuation in the xy plane, because the swinging fluctu- 
ation gives rise to {Sz{t)Sz{0)) . However, it is difficult 
to identify what kinds of magnetic fluctuations are asso- 
ciated with the in-plane and out-of-plane spin dynamics 
at the Ga site. To determine the origin of the tempera- 



ture dependence in G±{uj) and Gz{0), a theoretical study 
of the magnetic fluctuations with a 60° spin structure is 
desired. 

One of the characteristic features of the 1/Ti results 
on NiGa2S4 is that 1/Ti retains a large value down to 1 
K, i.e., far below Tf, which is also indicated by the /iSR 
relaxation rate in Fig. 10. The large value of 1/Ti is 
related to the fact that the NQR spectrum affected by the 
inhomogeneous static field cannot be observed between 2 
K and Tf. In an ordinary magnetic transition and spin 
glass ordering, 1/Ti shows a divergence at the magnetic 
ordering temperature Tm and decreases abruptly below 
Tm, and then the internal field emerges. Gontrastingly, 
in NiGa2S4, l/Ti is too short to be measured by the 
Ga-NQR measurements down to 2 K (~ 7f/5). This 
indicates that the magnetic fluctuations are not quenched 
abruptly but gradually decrease or stay constant down to 
2 K. 

Another characteristic feature of the 1/Ti behavior is 
that the enhancement of 1/Ti does not start just above 
Tf but at 80 K, which is far above Tf. Such behavior of 
1/Ti suggests that the release of the magnetic entropy 
at Tf is small. It is considered that the gradual decrease 
of the magnetic entropy is the reason why the anomaly 
of the specific heat at Tf is so broad that the anomaly 
seems to be different from ordinary magnetic ordering. 
We point out that the small change of the entropy at 
the magnetic anomaly is one of the characteristic fea- 
tures of frustrated magnetism, because a similar broad 
maximum in the specific-heat results is often observed in 
not only in triangular compound a^^i^° but also in kagome 
compounds.-^i25 

Furthermore, there remains a question to be under- 
stood when we interpret the temperature dependence of 
1/Ti on the basis of the spin- wave model. As reported 
from the neutron experiment, the magnetic correlation 
length in NiGa2S4 is short with only 6 lattice constants 
even at 1.5 K, which makes NiGa2S4 resemble more con- 
ventional spin glasses than the long-range ordered 2-D 
AF magnet. The presence of the spin wave seems to be 
unrealistic in such a short-range ordered compound. In 
addition, the spin-wave model is considered to have a 
significant field dependence of 1/Ti, which seems to be 
contradict with the field insensitive C(T).^ 

However, we point out that the similar situation was 
observed in SGGO,^'' in which the magnetic correlation 
length is only two times of inter-Cr spacing^^ but specific 
heat C(T) varies as T^ at low temperatures<^ In addi- 
tion, the similar slowing down of the Gr spin fluctuations 
when approaching to the susceptibility-cusp temperature 
Tg ~ 3.5 K was reported from the /iSR and Ga-NMR 
measurements.'*'^^ It is considered that the linear disper- 
sion indicated from the specific-heat behavior does not 
imply the presence of the long correlation length. We 
suggest the possibility of a cluster glass, which might be 
consistent with all existing results except for the field in- 
sensitive behavior of C(T). The spin- wave excitations 
arise inside the cluster and the clusters freeze out inde- 
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pendently. It is considered that the size of the cluster 
might be approximately 6 lattice constants. The per- 
sisting spin fluctuations below T{ might be associated 
with the spins at the cluster boundaries. To examine this 
scenario, ac susceptibility measurements and theoretical 
studies of the field dependence of N{uj) in the frustrated 
systems are highly desired. 

Quite recently, Tsunetsugu and Arikawa pointed out 
that spin nematic order can explain the experimental 
facts of absence of the magnetic long-range order, a 
power-law behavior C ^ T^, and the incommensurate 
wave vector q ~ (1/6, 1/6, 0).— The order parameters of 
the spin nematic order are not related to ordinary static 
spin dipole moments, but related to anisotropy of spin 
fluctuations. The present results of l/Ti ~ and the 
anisotropic spin fluctuations below T{ seems to be consis- 
tent with the spin nematic order. However, it is a crucial 
point whether the spin freezing behavior accompanied by 
the internal magnetic flelds below 2 K, observed by the 
NMR/NQR and /xSR measurements, is consistent with 
their scenario. 

On the other hand, Kawamura and Yamamoto studied 
the ordering of the classical Heisenberg antiferromagnet 
on the triangular structure with bilinear and biquadratic 
interactions.— They suggested that a topological phase 
transition at a flnite temperature driven by topological 
vortices although the spin correlation length remains fi- 
nite even below the transition point. They pointed out 
that the magnetic anomaly at T{ in NiGa2S4 might origi- 
nate from a vortex-induced topological transition related 
to Z2 vortices. In their vortex scenario, the magnetic- 
correlation time does not truly diverge at Tf, but only 
grows sharply at T{ exceeding the experimental time 
scale, and remains short over wide temperature range 
below Tf. Correspondingly, magnetic-correlation length 
^ is short and remains finite in this temperature range, 
and the onset of the magnetic long-range order, which is 
characterized by the exponential divergence of ^, is ob- 
served only at still lower temperatures The magnetic 
properties revealed with the present NQR, /iSR and the 
neutron experiments above 2 K seem to be consistent 
with the above vortex transition. In addition, if the tem- 
perature for long-range order in the vortex scenario is 
considered to be ~ 2 K, the observation of spin- wave-like 
excitations similar to those in a long-range ordered 2-D 
AF magnet can be understood. To examine this vortex 
scenario, the neutron experiments below 2 K are highly 
desired. 



V. CONCLUSION 

We have performed ss^^iQ^^.p^^j^/jv^Qj, ^gj, 
surements on triangular antiferromagnet NiGa2S4. NMR 
and NQR spectra above Tf indicate the existence of two 
Ga sites with different local symmetries, although there 
exists only one crystallographic site in perfect NiGa2S4. 
The intensity ratio of two NQR peaks is approximately 



4 : 1 and the linewidth of the intense peak is narrower 
than the other. At present, the origin of the two NQR 
peaks is not fully understood, but it is speculated that a 
tiny amount of sulfur disorder and/or planar defects such 
as stacking faults might induce the different Ga site. 

1/Ti and I/T2 were measured at both Ga sites and 
divergent behavior of 1/Ti and l/r2 was observed at Tf 
at both Ga sites. It was found that 1/Ti and I/T2 show 
the same temperature dependence above Tf , indicative of 
isotropic spin dynamics in the rapid-motion limit. The 
observed isotropic spin dynamics is considered to be char- 
acteristic of a Heisenberg spin system. With decreasing 
temperature, magnetic correlations start to develop be- 
low |6'w| 80 K, which is far above Tf . Spin fluctuations 
continue to slow down below and become nearly 

static below Tf = 10 K, where the specific heat C(T) 
shows a broad maximum. The remarkable divergence of 
1/Ti and I/T2, which is observed over wide temperature 
range between 0w and Tf , is characteristic of frustrated 
systems, because geometrical frustration suppresses the 
magnetic anomaly down to low temperatures. 

However, the wide temperature region between Tf and 
2 K, where the NQR signal was not observed, suggests 
that the Ni spins do not freeze immediately below Tf, 
but keep fluctuating down to 2 K with the MHz fre- 
quency range. Below 0.5 K, all components of 1/Ti follow 
a behavior. Below Tf, we found a broad spectrum 
and an inhomogeneous distribution of Ti in NMR and 
NQR measurements. These results indicate the freez- 
ing of magnetic moments and the emergence of inho- 
mogeneous static magnetism below 2 K. These are also 
suggested from the /iSR experiments. The relaxation 
rates decrease below Tf, and 1/Ti follows a depen- 
dence below 0.8 K, and the overall temperature depen- 
dence of 1 /Ti is roughly interpreted by the 2-D spin- wave 
model. In addition, on the basis of the observed differ- 
ence in the temperature dependence of 1/Ti and I/T2, 
the spin dynamics is interpreted to become anisotropic 
below 0.8 K. These experimental results strongly sug- 
gest that short-range magnetic order with incommensu- 
rability and/or inhomogeneous static moments is realized 
below Tf, which is consistent with the neutron experi- 
ments. We suggest that the w-linear dispersion implied 
from the speciflc-heat measurements is understood as a 
consequence of the inhomogeneous magnetic flelds due 
to the inhomogeneous short-range magnetic order. The 
spin dynamics in the short-range ordered state, which is 
considered to be fleld insensitive, and the magnetic cor- 
relations developing below 80 K are interesting issues in 
frustrated magnetism to be further studied from theoret- 
ical and experimental points of view. 
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